Generalized van der Waals theory for the twist elastic modulus and helical pitch of 

cholesterics 



H. H. Wensinlfl and G. Jackson 
Department of Chemical Engineering, Imperial College London, 
South Kensington Campus, London SW7 2AZ, United Kingdom 
(Dated: May 26, 2009) 

We present a generalized van der Waals theory for a lyotropic cholesteric system of chiral sphero- 
cylinders based on the classical Onsager theory for hard anisometric bodies. The rods consist of a 
hard spherocylindrical backbone surrounded with a square-well potential to account for attractive 
(or soft repulsive) interactions. Long-ranged chiral interactions are described by means of a simple 
pseudo-scalar potential which is appropriate for weak chiral forces of a predominant electrostatic 
origin. Based on the formalism proposed by Straley [Phys. Rev. A 14, 1835 (1976)] we derive 
explicit algebraic expressions for the twist elastic modulus and the cholesteric pitch for rods as a 
function of density and temperature. The pitch varies non-monotonically with density, with a sharp 
decrease at low packing fractions and a marked increase at higher packing fractions. A similar 
trend is found for the temperature dependence. The unwinding of the helical pitch at high densities 
(or low temperatures) originates from a marked increase in the local nematic order and a steep 
increase of the twist elastic resistance associated with near-parallel local rod configurations. This 
contrasts with the commonly held view that the increase in pitch with decreasing temperature as 
often observed in cholesterics is due to layer formation resulting from pre-smectic fluctuations. The 
increase in pitch with increasing temperature is consistent with an entropic unwinding as the chiral 
interaction becomes less and less significant than the thermal energy. The variation of the pitch 
with density, temperature and contour length is in qualitative agreement with recent experimental 
results on colloidal fd rods. 

PACS numbers: 83.80.Xz, 61.30.Cz, 82.70.Dd 



I. INTRODUCTION 

Cholesteric liquid crystals (LCs) display a variety of ex- 
traordinary features due to the existence of a mesoscopic 
helical structure which is best known for its exception- 
ally large optical rotational power employed in LC dis- 
play technology. In contrast to a common nematic phase, 
where the nematic director is homogeneous throughout 
the system, the cholesteric (chiral nematic) phase is char- 
acterized by a helical arrangement of the director field 
along a common pitch axis. As a result, the cholesteric 
phase possesses an additional mesoscopic length scale, 
commonly referred to as the 'pitch length', which char- 
acterizes the distance along the pitch axis over which the 
local director makes a full revolution^. 

Derivatives of cholesterol, chiral molecules which were 
the first mesogenic substances to be recognized in stud- 
ies of the melting point and optical properties of car- 
rot extracts by Reinitzer— and LehmannSi, belong to the 
thermotropic class of liquid crystals where phase transi- 
tions are brought about by a variation in the temper- 
ature; a large number of thermotropic mesogens with 
a low to moderate molecular weight that form chiral 
nematic phases have now been isolated or synthesized. 
Their widespread use in optoelectronic applications (e.g., 
twisted nematic liquid crystal displays in laptop com- 
puters, televisions, and mobile phones) is a direct con- 
sequence of the unique rheological, electrical and optical 
properties imparted by the chiral structures. 

Lyotropic chiral systems, involving high molecular- 



weight particles in solution where the ordering behav- 
ior is primarily governed by the solute concentration, 
are also common. Examples are bio-colloidal systems 
such as DNAii^ and the rod-like /d-virus^, stiff poly- 
mers such as polypeptidesii^, polysaccharides^ or cellu- 
lose derivatives^^, and chiral micellesii. In these sys- 
tems, the cholesteric pitch is very sensitive to con- 
centration, temperature as well as the solvent condi- 
tions such as ionic strength and pH. The dependence 
of the pitch upon these variables has been the sub- 
ject of intense experimental research (see for example 

^^fj,^^j^j,^,^4 .12.13.14.15.16.17.18.19.20.21.22 -| 

Theoretical attempts to predict the behavior of the 
cholesteric pitch are challenging owing to the complexity 
of the underlying chiral interaction^^ and the inhomo- 
gencous and anisotropic nature of the phase. Course- 
grained model potentials aimed at capturing the essen- 
tials of the complex molecular nature of the electrostat- 
ics of the surface of such macromolecules have been de- 
vised mainly for DNA^ii^iS^i^iS^. A more general elec- 
trostatic model potential for chiral interactions was pro- 
posed much earlier by Goossens^S based on a spatial ar- 
rangement of dipole-dipolc and dipole-quadrupole inter- 
actions which can be cast into a multipole expansion in 
terms of tractable pseudo-scalar potentials^. This type 
of electrostatic description of the chiral interaction can 
be combined with a Maier-Saupe mean-field treatment 
(see, for example reference a'^^i'^^'^'^i'^^'^^i'^^ , or with a bare 
hard-core model and treated within the seminal theory of 
Onsager— (as in the recent study of chiral hard sphero- 
cylindcrs for systems with perfect local nematic order~) . 
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Taking the alternative view of a steric origin for the 
chiral interactions, Stralej*^^ combined a leading order 
pseudo-scalar form with a hard rod model to create 
a basic chiral model (hard threaded rod) for lyotropic 
cholesterics. Though influenced by the work of Goossens 
on electrostatic forces, the potential used by Straley 
is appropriate for short-ranged steric chiral interactions 
mediated by a thin helical thread enveloping each rod. 
By extending the theory of Onsager— to systems with 
non-uniform director fields, microscopic expressions for 
the macroscopic twist energy associated with a bulk 
cholesteric structure can be deduced. The theoretical 
treatment was later elaborated in a paper by Odijk^ 
leading to microscopic scaling expressions for the twist 
elastic constant and helical pitch of both rigid and semi- 
flexible rods. Similar relations, albeit with different scal- 
ing exponents, were obtained by Pelcovits^ based on a 
corkscrew model. 

At this stage one should also acknowledge the related 
studies on the link between the orientational order pa- 
rameters and the elastic constants in nematics, which 
have been influential in shaping some of the theories 
of chiral phases. A mean-field treatment of particles 
with electrostatic chirality but no shape anisotropy has 
been developed^>^i^, as have Onsager DFT theories for 
anisotropic hard rod-like particles^i^i^>^ and hard rods 
with attractive mean-field a^^'^° , the latter being closely 
related to the achiral contribution of the model employed 
in our study. The principal finding of this body of work 
is that the elastic constants are predicted to be propor- 
tional to the square of the nematic order parameter, at 
least for system with weak orientational order, confirm- 
ing experimental observation^ and the results of molec- 
ular simulation^. Such a treatment has also been ex- 
tended to a description of the elastic constants in smectic 
phases^^^. 

The underlying microscopic physical feature (steric or 
electrostatic) responsible for the formation of chiral ne- 
matic phases is still a matter of debate and controversy. 
In has been known for some time that the nature of 
the solvent can have a dramatic effect on the pitch of 
cholesteric phases, and can even reverse the sense of 
the twist (e.g., see the study of Robinson^ on poly-7-L- 
benzyl-L-glutamate PBLG, a synthetic polypeptide with 
an a helical conformation, in achiral solvents such as 
dioxane and dichloromethane), hinting to a solvent ef- 
fect of electrostatic origin. In a beautifully revealing but 
rather overlooked paper, Coates and Gray^ showed that 
the replacement of a hydrogen by a deuterium atom on 
the carbon backbone of an originally achiral thermotropic 
mesogen is sufficient to induce a cholesteric structure; the 
carbon-hydrogen and carbon-deuterium bond lengths are 
both 1.085 A so one would expect the "steric shape" of 
both molecules to be very similar, indicating that in this 
case at least the chirality is of a weak and subtle electro- 
static nature. 

The fact that the pitch of the cholesteric phase in aque- 
ous solutions of the filamentous fd virus is very sensitive 



to the ionic strength of the medium but still persists af- 
ter coating the virus with polyethylene oxide polymer 
(which would mask any short-range chirality in the par- 
ticle shape) is also indicative of an electrostatic origin 
to the interparticle chiral interaction in such lyotropic 
system a^^i^^ . This lends credence to the use of electro- 
static interactions of the type proposed by Goossens^^ 
decorated with an achiral non-spherical core as physi- 
cally reasonable microscopic models for chiral systems; 
an additional repulsive steric chiral core could of course 
also be incorporated in a more realistic treatment (e.g., 
see reference^), but this is beyond the scope of our work. 

The helical pitch of cholesteric thermotropic meso- 
gens is almost invariably found to be a decreasing func- 
tion of temperature^, which is commonly attributed to 
the presence of an underlying smectic-A phase at lower 
temperatures: a twisted chiral structure is incommen- 
surate with the smectic layering leading to an unwind- 
ing of the pitch as one approaches the transition^^i^. 
What is surprising is that this thermally induced de- 
crease in the pitch can occur over many decades in tem- 
perature, where one would not expect effects due to pre- 
transitional smectic order. An increase in pitch with in- 
creasing temperature has been reported for the choles- 
terol ester, cholcsteryl [2-(2-ethoxyethoxy-ethyl] carbon- 
ate (CEEC)^; interestingly in the case of CEEC the 
transition from the isotropic phase to the chiral nematic 
does not appears to be followed by a transition to a smec- 
tic phase with a further decrease in temperature, but 
rather to the formation of a crystalline stated. A re- 
markable sense inversion in the helical pitch with temper- 
ature has also been found in thermotropic (solvent free) 
polypeptides^ and cellulose derivatives^, and in mix- 
tures of right-handed cholesterol chloride and left-handed 
cholesterol myristate^, indicating a subtle balance in the 
forces giving rise to chiral phases. 

The situation is just as intriguing in the case of ly- 
otropic systems where depending on the range of temper- 
ature, both negative and positive slopes of the pitch can 
be observed. At relatively high temperatures a marked 
increase of the pitch has been found in polypeptide sys- 
tems (mixtures of PBLG in dioxane, chloroform, and 
dichloromethane)^ and in aqueous solutions of /d-virus 
rod&i^. In the latter system the pitch is found to increase 
with increasing concentration well before the transition 
to a smectic phase. An unwinding of the cholesteric 
phase (with a corresponding increase in the pitch) is 
also observed in aqueous solutions of DNAi^ as the con- 
centration of the macromolecules is increased towards a 
high-density hexagonal columnar (positionally ordered) 
stated; this is analogous to the divergence of the pitch 
found in thermotropic mesogens as one approaches the 
smectic state on lowering the temperature. 

The early molecular-field approaches based on the 
Maier-Saupe electrostatic picture of the mesogenic 
interactionSii^ fail to give a consistent picture of the 
variation of the pitch with temperature for cholesteric 
phases, mainly because of the lack of a hard-core exclude 
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volume contribution in the interaction which plays an 
essential role in the stabilization of the orientationally 
ordered phased. 

The important role of the mesogen's shape in under- 
standing the temperature dependence of the pitch was 
pointed out early on by Kimura et al^. The insensi- 
tivity of the pitch to temperature found in lattice sim- 
ulations of sites interacting through a Goossens electro- 
static potential^ supports the view that one requires a 
balance between the opposing twist elastic forces (which 
are primarily a consequence of the repulsive interactions) 
and the chiral torque (due to the cholesteric interactions) 
to observe the subtle dependencies found for the pitch. 
A proper account of the rod-like backbone turns out to 
be essential to explain the unwinding of the cholesteric 
structure with increasing temperature^^i^. 

It is clear from the preceding discussion that a num- 
ber of important questions remain unresolved in our un- 
derstanding of the origin of chirality and the related 
dependence of the helical pitch on the thermodynamic 
variables such as temperature and composition. Docs 
one require the presence of an underlying smectic (or 
columnar) phase to observe an increase in pitch with de- 
creasing temperature? Under what conditions would one 
expect a thermally induced unwinding of the chiral nc- 
matic phase and when would a non-monotonic temper- 
ature dependence of the pitch be obtained? What are 
the competing roles of repulsive excluded volume inter- 
actions, and isotropic and anisotropic (chiral and achiral) 
attractive interactions in stabilizing a twisted equilibrium 
structure? As de Gennes and Prost state in their mono- 
graph when referring to the temperature dependence of 
the pitch: "Their origin is not yet quite clear" . In this 
paper, we will build upon the idea of coupling a hard core 
with chiral electrostatic interactions and present a micro- 
scopic theory for rigid chiral spherocylinders of arbitrary 
aspect ratio. 

Our motivation for the analysis is three-fold. First of 
all, most investigations thus far have been restricted to 
infinitely thin rod-like species. We show that the rod 
thickness plays a crucial role in the behavior of the pitch 
in cholesteric systems with strong orientational order and 
helps to account for a non-monotonic dependence of the 
pitch with temperature and concentration. Secondly, 
the role of long-range chiral and achiral dispersive forces 
are investigated here by means of a simple square-well 
(SW) potential. This allows us to probe the generic tem- 
perature dependence of the pitch for chiral rods where 
additional achiral attractions (or soft-repulsions) are at 
play. Finally, we pay particular attention to the impli- 
cations of chirality in the intcr-particlc interactions on 
the phase behavior. Our results are entirely algebraic 
and contain explicit expressions for the twist energy, elas- 
tic modulus, and cholesteric pitch as a function of den- 
sity and temperature for a given rod aspect-ratio and set 
of intcrmolccular SW parameters such as the interaction 
range and the relative magnitude of the chiral and achi- 
ral attractive/soft-repulsive rod interactions. The results 



are fully consistent with experimental results for /d- virus 
rods and may prove helpful in interpreting observations 
in other cholesteric systems. 

This paper is laid out as follows. We start with a brief 
exposition of Straley's theory for the deformation free en- 
ergy of the cholesteric state in Sec. II. A suitable sphe- 
rocylindcr potential is introduced in Sec. Ill and incor- 
porated into the deformation free energy through a van 
der Waals treatment based on the Onsager-Parsons the- 
ory for anisometric hard bodies. In Sec. IV, microscopic 
expressions for the twist parameters and cholesteric pitch 
is derived. These are analyzed and the predictions com- 
pared with experimental results in Sec V. Finally, some 
concluding remarks are made in Sec. VI. Technical issues 
will be relegated to a number of Appendices. 

II. MICROSCOPIC THEORY FOR THE 
TWISTED NEMATIC 

To describe the properties of the cholesteric phase we 
will follow closely the analysis proposed by Stralcy22i^. 
The aim is to calculate the distortion free energy as- 
sociated with non-uniform nematic director fields. In 
a distorted (e.g., twisted) nematic phase the director 
is no longer spatially uniform but depends on position. 
The same holds for the orientational distribution func- 
tion (ODF), /(u • n(r)) which describes the probability 
of finding a particle with a given orientational unit vector 
u with respect to a locally varying director n(r). Within 
an Onsager-type formulation^, the excess Helmholtz free 
energy of a distorted nematic state of N particles in a 
volume V can be cast in the following form: 

= ^ / ^s*(ri2; ill, fi2)/(ai • fi(ri))/(u2 • ft(r2)), 

(1) 

with J ds ~ JJJJ (iri(ir2(iuidu2 and p ~ N/V the num- 
ber density. The kernel $ accounts for the pair interac- 
tion of two rods for a relative centre-of-mass separation 
ri2 ~ r2 — ri and orientations Ui, U2. This quantity 
will be fully specified later. Note that, strictly, / — 1/Att 
in the isotropic state where the particle orientations arc 
completely random and the director field becomes irrel- 
evant. 

If the spatial variation of the local nematic director is 
weak such that the associated distortion wavelength is 
much larger than the particle dimensions, the ODF may 
be approximated by a Taylor expansion. For the spatial 
integration it is expedient to switch to a new coordinate 
system, R±ri2/2 for i = 1,2, in terms of R = (ri -|- 
r2)/2 and the centre-of mass distance ri2. Expanding the 
ODF then gives 

/(u..n(r,)) = /(u, • fi(R)) ± ^i^/(u, . n(R)) 

+0 [(V • rr2f] ■ (2) 
With this result, the excess free energy for weak director 



4 



gradients becomes after some rearrangement 

P"^ = ^ j ds'<I>(ri2;ui,U2)/(iii-n(R))/(u2-n(R)) 

+ ^ y'ds'<I>(ri2;ui,U2)/(ui.n(R))/(u2-n(R)) 

x[(ri2 • V)n(R) -112] 
- ^ /ds'<I>(ri2;ui,U2)/(ui.n(R))/(u2-n(R)) 



[(ri2 • V)n(R) • ui][(ri2 • V)n(R) • U2] 



(3) 



with / ds' = //// (iRdri2(iui(iu2 and / = df{ui ■ 
n(R))/9(ui-n(R)). Let us now consider a twist deforma- 
tion of the director with a heUx axis afong the z-direction 
of the laboratory frame. Assuming the focal nematic di- 
rector to describe a perfect helix, we can parametrize the 
director field as follows 



n(R) = cos{qZ)Sc + sm{qZ)y, 



(4) 



with q = 2t:/p the magnitude of the pitch wave vector 
and p the length of the cholesteric pitch. For small wave 
vectors, ga <C 1 (with a the typical range of the pair po- 
tential) the trigonometric functions can be expanded to 
leading order and the spatial dependence of the director 
approximated as 



n(R)=x + gZy + 0(<j2). 



(5) 



The free energy density of the twisted nematic state can 
then be expressed as 



^duidu2A/o(ui,U2)/(ui)/(u2) 



-Ktq+-K2q\ 



(6) 



where the coefficients pertain to the twist energy and 
twist clastic contributions, respectively: 



A't(n-Vxn) = jj duidu2Mi{ui,U2) 

X-U2y/(U1)/(U2) 



K2{n • V X n)2 



rfUidU2A/2(ui, U2) 



XUiyU2y/(ui)/(u2). 



(7) 



(8) 



The first contribution, (Eq. ([71)); is non- vanishing only if 
the rod interactions are chiral, as we will see later on. 
It provides a measure for the torque-field exerted by the 
microscopic chiral interaction which favors the twist dis- 
tortion. The second term, Eq. ([8]), provides a microscopic 
expression for the twist elastic modulus K2 of the Frank 
elastic free energy^. The elastic contribution is a restor- 
ing (Hookian) term in the free energy which resists the 



twist deformation. The quantities Mi(ui,U2) are given 
by the following moment spatial integrals over the inter- 
action kernel <J>: 

A4(ui,U2)= y"dri2$(ri2;ui,U2)zi^2, ik = 0,l,2). 

(9) 

In case of achiral hard particles, only the reference term 
in Eq. ([6|) needs to be considered. Moreover, with 
A/o(ui,U2) = Wcxci(ui, U2), the excluded volume of a pair 
of particles, the original Onsager free energy^ is recov- 
ered as required. 

For chiral nematic systems the equilibrium value of the 
cholesteric pitch is found by balancing the chiral forces in- 
ducing the twist deformation with the elastic forces which 
favor the nematic. Minimizing the total free energy with 
respect to the local ODF and q leads to following coupled 
set of stationarity conditions: 



F 



id 



V 



- duf{u)^ = (10) 

(11) 



d fF 



F" 



V 



= 0, 



where /i is a Lagrange multiplier associated with the nor- 
malization constraint for the ODF (J o?u/(u) = 1) and 



k^T J dup/(u)[lnVp/(u) - 1] (12) 



is the exact ideal free energy of a spatially uniform sys- 
tem (V is the thermal volume of the rod). For weak twist 
deformations it is safe to assume that the local ODF re- 
mains unaffected by the twist and that it adopts the same 
form as in the nematic. Denoting /(u) = /o(u), the equi- 
librium ODF of the nematic phase, the equilibrium value 
for the pitch wave vector at a given density p is given by 
the ratio of the average chiral and elastic forces^: 



q = 



Ktifo] 
K2U0] 



(13) 



The expressions for the twist parameters can be made 
analytically tractable by using a Gaussian trial function 
Ansatz to describe the local ODF— This is done in 
Sec. IV. First we have to specify the form of the inter- 
action potential and kernel $ based on a suitable pair 
potential for chiral spherocylinders. 



III. GENERALIZED VAN DER WAALS 
THEORY FOR SW RODS 

Whilst the seminal view of Onsager— that the repul- 
sive inflexible core of a particle gives rise to orientation- 
ally ordered phases is now very well established, the spe- 
cific nature of the dispersive and polar interactions can 
have an important influence on the macroscopic struc- 
tures that are observed. For example, in the case of hard 
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rods with central point dipoles, layered liquid crystalline 
phases such as the sniectic-A are favored, while the effect 
on the isotropic-nematic transition appears to be small 
and in some cases unfavorabl o^^i^^'^'^ ; for molecules with 
terminal dipoles, the nematic phase is stabilized relative 
the smectic phase^i. As can be inferred from the discus- 
sion in the introductory section, the influence of attrac- 
tive interactions is all the more beguiling and subtle in 
the case of systems with chiral interactions, where the 
pitch of the helix is found to be very sensitive to the 
balance of forces that the molecules experience. 

Let us consider an ensemble of hard spherocylinders 
(HSC), cylinders of length L capped by hemispheres of 
diameter D. The rods can be rendered chiral by intro- 
ducing the simple pseudo-scalar potential, proposed by 
Goosscna^. For any non-overlapping configuration of a 
pair of rods the chiral contribution reads 



chiral 



-£212M(ri2)r212(fl2; Ui, U2), 



(14) 



which consists of a radial part u{r 12) and an orientation- 
dependent pseudo-scalar defined as 



T212(fl2;Ui,U2) = (Ul • U2)(U1 X U2 • V12), 



(15) 



which in fact represents the first non-trivial term in a 
series expansion in terms of generalized chiral pseudo- 
scalars T2i{2k-i)2j ■ For the present case of weak chiral 
interactions it suffices to retain only the first term. As 
T212 changes sign upon interchanging particle positions, 
ri2 — fi2, while keeping the orientations fixed, the 
pseudo-scalar imparts a chiral interaction. The sign of 
the amplitude £212 defines the handedness of the chiral 
interaction and the corresponding helical mesostructure. 

In Goossens' model^^ the radial part decays steeply 
via u{r 12) = ^/ 1^121 which arises from a summation 
over electrostatic dipolar interaction sites located on each 
rod. Here, we employ a much simpler dependence based 
on a simple SW form with range A. We thus specify 
u(ri2) = H{X — 7'i2), with H a Heaviside step function. 
To take into account the effect of achiral attractive (or 
soft-repulsive) forces we introduce an additional achiral 
SW potential with amplitude (well depth) eooo • For sim- 
plicity, we assume both the chiral and achiral SW poten- 
tials to have the same interaction range A: 

achiral = — e000^^(A — ri2) 
$chiral = -e2127212(fl2;Ui,U2)i?(A - ri2), (16) 

for any non-overlapping rod pair configuration. Putting 
all contributions together, we arrive at the following total 
pair potential for chiral SW rods: 



oo 



^tot (ri2;ui,U2) 



"^000 







ri2 < a 

J212?212 cr < 7'i2 < A 

ri2 > A, 

(17) 

with o'(fi2; Ul, U2) the centre-of-mass contact distance 
between two hard spherocylinders at given (relative) ori- 
entations. Henceforth, we shall fix the SW range at 



X = L + D (the so-called "square peg in a round hole" 
model). It is advantageous to introduce a reduced tem- 
perature defined as T* ~ fcsT/leoool- An expression for 
the free energy of the nematic phase of particles consist- 
ing of a hard anisotropic core with attractive interactions 
can be obtained from a first-order perturbation theory 
around a suitable hard-core reference free energy. A gen- 
eralized van der Waals (GvdW) form can be expressed 



a* 



75.76.77 



-^GvdW 

V 



^k^TG{^) j I rfUidU2/(Ui)/(U2) 



12 / dri2rl2 



X I dr 

- yJJ dnidn2fiui)fiu2) 

X J dri2 J <iri2ri2 (£000 +£2127212). (18) 

The first contribution is the Onsager-Parsons excess 
free energy accounting for the hard-core repulsive part 
of the pair potential. It is based upon a scaled sec- 
ond virial approximation according to the Parsons-Lee 
approach22iZ2ii22. It involves a mapping of the radial dis- 
tribution function for the anisotropic particles onto that 
of an equivalent hard-sphere system via the virial equa- 
tion. The rod free energy can ultimately be linked to 
the Carnahan-Starling^iiS^ expression for hard spheres 
which provides a simple strategy to account for the ef- 
fect of higher-body interactions, albeit in an implicit and 
approximate manner. The scaling factor 



G(</>) = 



1 



(1-, 



(19) 



reduces to unity in the Onsager limit L/D — > oo where 
the packing fraction of the nematic phase at the order- 
ing transition vanishes. The second term of Eq. (|18p is 
the contribution to the free energy due to the attractive 
interactions, at the mean-field level of description (i.e., 
any correlations in the particle positions are neglected). 
This type of augmented van der Waals equation of state 
is commonly employed in studies of homogeneous fluids 
(e.g., see references^i^) , and can even be used to de- 
scribe complex fluid phase equilibria in a quantitative 
manner (for instance the liquid-liquid phase behavior of 
hydrocarbons and perfiuoroalkanes^, or aqueous mix- 
tures of hydrocarbons^^ and amphiphiles^i) . 

If the GvdW excess free energy is mapped onto the 
general form given by Eq. ([T]) we can obtain the following 
expression for the total interaction kernel 



$(ri2;ui,U2) 



kBTGi^) ri2 < a 

— £000 ~ £2127212 C < 7*12 < A. 



.... (2°) 

It is important to note that $ is not invariant with re- 
spect to ri2 — * — ri2 as would be the case for achiral 
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FIG. 1: Phase diagram of achiral (£212 = 0), attractive SW spherocylinders (eooo > 0) of aspect ratio L/D = 10 and range 
A — L+D. The phase behavior is represented in (a) the reduced temperature T* — feBT/|eooo| versus molecular packing fraction 
4) = vop plane with vo = {'k/A)LD^ + {tv/6)D^ the spherocylinder volume, and (b) the reduced pressure P* = Pvo/\£ooo\ versus 
T* plane; the continuous curve in (a) represents an isobar for P* — 120. The calculations are based on the Onsager and 
Gaussian trial functions, Eq. (|24|l and Eq. (|25|l . respectively, (c) Variational parameter a of the coexisting nematic phase 
versus T* . 



interactions. In the case of ranges of the SW interac- 
tion that are at least as long as the longest dimension of 
the particles (L + D) the integration of the interaction 
decouples neatly into a repulsive and attractive contribu- 
tions; the integration is not as straightforward for shorter 
ranged interactions. 

The fluid phase behavior for achiral attractive rods 
(^212 = 0) has been discussed extensively in Rcfs. [tgUttI 
and a typical example is reproduced in Fig. [T] For this 
calculation we only need the zeroth moment Mq which 
readily follows from Eq. ^ and Eq. ([^ : 

A/o(ui,U2) = (fcBTG(</.)+£ooo)wcxcl(7)-£oooYA^ (21) 

where Woxci is the excluded volume between two hard 
spherocylinders at a relative angle 7 = arcsin |ui x U2I: 



Woxcl(7) 



120" 



47r 



2L^D\au-i-f\ + 27rD^L+ —D^. (22) 

o 



The ODF /Ar(u) of the nematic state is formally given 
by the solution of the self-consistency equation emerging 
from Eq. HI 



fN in) = 



cxp [~pJdu'Mo{u,u')fN{u')] 
J ducxp [-P J du'A/o(u, u')fN{n')] 



(23) 



which is not amenable to further analysis and is solved 
numerically. The double orientational averages appear- 
ing in Eq. can be made analytically tractable by 
adopting a simple algebraic trial form for the ODF. In his 
original paper Onsager— introduced the following form 
to describe the distribution of angles in the uniaxial ne- 
matic state: 



fo{0) 



acosh(acos( 
47r sinh a 



(24) 



where cos = u • n is the polar angle between the axis 
of symmetry of the particle and the director and a > 



a variational parameter measuring the degree of nematic 
order. Note that a = in the isotropic phase. The 
implications of Eq. on the GvdW free energy have 
been studied in detail in Ref . [ttI . A simpler trial function 
has been proposed by Odijk and Lekkerkerker— , based 
on a Gaussian: 



fG{0) 



in 



exp 



0<e<TT/2, (25) 



and its mirrored version for the interval tt/2 < 9 < n. 
Unlike Eq. , the Gaussian trial form does not reduce 
to the correct isotropic constant l/47r for a = and we 
must therefore require a 3> 1 for reasons of consistency. 
In fact, Eq. ([24]) becomes identical to the Gaussian form 
in the asymptotic limit a ^ 00. The results from both 
distributions will therefore be virtually indistinguishable 
for large a, as is illustrated in Fig. \Vp. 

Let us quote the following Gaussian averages^S: 

(ln/(u)) - ln47ra- 1 

((sin 7)) ~ (lY^' (a»l), (26) 

where the brackets denote orientational averages: 



dufaiu), 



((•)) ^ JJ dUidU2/G(Gl)/G(G2). (27) 

The nematic order parameter S can be approximated 
by: 



S= (P2(cos0)) ~ 1- -, 
a 



(28) 



with V2 the second-order Legendre polynomial. Using 
the asymptotic expressions from Eq. ((26|) in the ideal 
and excess free energy yields a simple algebraic expres- 
sion for the free energy. The corresponding phase equi- 
libria can be analyzed without difficulty and the resulting 
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phase diagrams are shown in Fig. [1] the coexisting den- 
sities are obtained by numericaUy solving the equality 
of pressure Pi — ~{dF/dV)N.T and chemical potential 
Hi = —{dF/dN)v,T of each phase i. At low to mod- 
erate temperatures, a coexistence between isotropic gas 
(Ig) and liquid (/l) phases is found in addition to the 
isotropic-nematic phase separation seen at higher den- 
sities. At the triple point temperature T* — 12.254 
the isotropic liquid binodal meets a triphasic Iq-Il-N 
equilibrium line. For systems with chiral interactions 
(^212 7^ 0) the phase behavior is altered by the twist con- 
tributions which we shall examine in the next Section. 



where ^ dri2 = / df 12 Jq dri2rf2 denotes a spatial 
integral over the spherocylinder excluded volume and 
dri2 — J dri2 rf?'i2^i2 a spatial integration over the 
SW range. All contributions of C(e2i2) are of negligible 
importance for weak chirality; more specifically, the ratio 
of the chiral and non-chiral interactions must be small, 
i.e.. 



^212 



£000 



< 1, 



(30) 



IV. CALCULATION OF THE TWIST 
PARAMETERS 

The central task in the description of chiral nematic 
phases within our GvdW theory is the calculation of the 
moment integrals pertaining to the twist energy and elas- 
tic modulus in Eq. ([7]) and Eq. ([8]) using the explicit inter- 
action kernel defined in Eq. (PO)) . Details of the specific 
treatment are provided in the next subsections. 



A. Twist elastic modulus 



where one should also note that G{4>) > 1. The chirality 
parameter Ec is, in principle, fixed by the detailed molec- 
ular structure of the rod-like mesogen and is expected to 
be small for most common chiral substances. In the case 
of strong chiral interactions, [cc ~ C'(l)], the twist elas- 
tic response will be affected by T212 (and higher order 
pseudo-scalar potentials) which severely complicates the 
analysis. The second integral in Eq. ([M)) is easily eval- 
uated by exploiting the symmetry of the SW potential 
along the z-axis and using cylindrical coordinates: 



The calculation of the twist elastic modulus consists 
of two steps. First, an explicit expression for the second 
moment M2 in Eq. ^ is required. Then, a double ori- 
entational average has to be carried out using an appro- 
priate form for the ODF according to Eq. ([5]). Inserting 
Eq. ((20|) into Eq. ([9]) and rearranging terms leads to 

M2(ui,U2) = iksTGicj}) + eooo) / c?ri2Z?2 

^ Ucxcl 

-eoon J dri2zj2+ 0(6212), (29) 



rfri2^i2 



27r 



dzi2zl2 



-A 



drr = 



47rA^ 
IE' 



(31) 



where r = (x^2 + Uvi) ■ The first term in Eq. 
involves a weighted spatial integral over the excluded- 
volume manifold spanned by two spherocylinders at fixed 
orientations. Details of this calculation can be found in 
Appendix A. The final expression for M2 reads: 



Af2(iii,£i2) = ikBTG{(b) 
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-eooo) [L^D\smj\ -{AI + AI + B^) 



LD' {—{A1C2- A2G1] 



8"' / x2 



Ai 



c? 



)} 



47r£oooA^ 
15 



(32) 



where cx 0{D^) is given by Eq. in Appendix A. 
The second moment excluded volume is independent of 
the chiral interaction in the limit of infinitesimally small 
twist distortions considered here. In Eq. ([32|) . Ai, B and 
Gi are orientationally dependent dot products specified 
in Eq. ((65)) of Appendix A. 

The elastic modulus is obtained from a double orienta- 
tional average of A/2 for which we shall invoke a Gaussian 



trial function for the ODF, c/. Eq. (|25|) . In principle a 
twist deformation of the director field breaks the uniaxial 
symmetry of the local ODF, which would now take on a 
biaxial form, and a suitable generalization of the Gaus- 
sian distribution involving an explicit dependence on the 
azimuthal angle would therefore be required. To keep 
the theory tractable, we shall ignore local biaxiality and 
use the original form of Eq. . As the degree of biaxial 
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nematic order is very small in the weak deformation limit 
considered here, it is unlikely to have a significant effect 
on the mesoscopic properties of the cholesteric phase. 

Since the Gaussian ODF is only appropriate if the ori- 
entational distribution is strongly peaked around the ne- 
matic director we may perform an asymptotic expansion 
of the Gaussian integrals to extract the leading order 
contributions for large a. Let us fix z = {0,0,1} and 
introduce 

Ui ~ {cos Qi , sin 6*,; cos i-pi , sin 6',; sin (p,; } , (33) 

the orientational imit vector of rod i = 1,2 in terms 
of the polar (di) and azimuthal (ipi) angle with respect 
to the nematic director which we have fixed along the 
x-direction of the Cartesian frame n = x = {1,0,0}. 
Expanding the dot products in Eq. ((65|) for 0^ <C 1 we 
obtain up to leading order: 

|sin7| - |7| - (6*2 + 6*2 -26*1 6*2 cos Av?)^/^ (34) 
A, - (i/2)6l,sinv3., 

B ^ _D(02 COS(y92 — 01 COS(y9i)/|7| 

C1-C2 ~ i^(0isin^i -02sin(p2)/|7l, (35) 

with A(p = kpi — Lp\. The twist elastic modulus can be 
recast into a finite series in terms of the inverse aspect 
ratio X = DjL: 



where 



K2D 

knT 



(36) 



with c = pL^D the dimensionless rod concentration and 



e = G(0)±^, 



(37) 



where (+) applies to an attractive square-well (eooo > 
0) and (— ) a soft-repulsive square-shoulder (eooo < 0) 
potential. The coefficients Hk represent double Gaussian 
averages involving the following angular quantities up to 
leading order in a 3> 1: 
2 

i^o ~ ^((e|7|(0?sinVi+^?isin2(/52))) + --- 



Hi 



sm if I 



: sm' 



>2))) 



H2 ~ y ((e(02COS(^2-0iCOS(^i)V|7|)) 

+ — {{Q{di sin - 62 sin^2) /I7I)) + 

2 

-ffa ~ ^{{g{02cosf>2-eicosipifh^)) 

+ — {{g{0i sin ^1 - 02 sin (^2)^7^)) + 



Hi 
H, 



15 

27ra2 

15 



((^^))=0, 



(38) 



Q ^6162 cos Lpi cos ip2 



(39) 



The brackets represent the following four-fold angular in- 
tegral: 



((•» 



n 



dOiOi exp 



27r 



(40) 

where we have used fa — 9/g/9(cos0) ~ olJg- It is 
evident that vanishes upon integration over the az- 
imuthal angles, irrespective of the form of the (uniaxial) 
ODF. A little inspection shows that the leading order 
contributions to Hi and are also eliminated by the 
double azimuthal integration. The remaining terms are, 
in principle, nonzero because I7I depends non-randomly 
on Aiy9. The results can be greatly simplified by changing 
to the new azimuthal variables ipi = and ip2 = "0 + Aiy9. 
The integration over ip can be carried out without diffi- 
culty. After some algebra the averages reduce to: 



Ho 
H2 

H. 



96 



{{eie2{el + 9lMcos/^v)) 



^((^1^2 [{91 +el)7 cos /^^ 
-(4+10cos2A</7)]/|7|)) + • 

— ((0102 C0SA(^)) + --- 



(41) 



Here, the brackets now denote Eq. ([40]) with the dou- 
ble azimuthal integral over ipi 2 replaced by a single one 
over the remaining angle A(^. Clearly, the azimuthal av- 
erage yields i?3 = so that we need only evaluate the 
two even contributions and H2. The first term cor- 
responding to infinitely long rods has been analyzed by 
Odijk in Ref. Issl . Eliminating cos A(/3 via the asymptotic 
expression for | sin7| [Eq. ([M]) ] leads to: 



Ho 



a" 
96 



m'ol)) + m6l{el + el)))] (42) 



We may now use the Gaussian averages in Appendix B 
to arrive at the compact expression. 



7 

"192 \a) 



(43) 



Similarly, by applying Eq. p4p to H2 we may simplify: 



H. 



a 
24 
1 

+ 8 



13 
24 



((i7i^?))-^((i7n) 



{{om-eDM)) 



(44) 



The first two Gaussian averages arc given in Appendix B. 
The last one cannot be calculated analytically, but the a 
dependence is easily established from a simple analysis of 
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the scaling, while the pro-factor can be found numerically 
(see Appendix B). The final result is: 



Ho 



(45) 



with K = 0.036926. For later reference, we also give 
the next leading order a-contributions to Hk- It can be 
shown that the odd terms Hi and H^ in Eq. ([55]) are fully 
eliminated (to all order in a) by the double azimuthal in- 
tegration. The higher order a-contributions to Hk (even 
k) are in principle non- vanishing and can be estimated by 
considering the angular quantities in Eq. (|38p taking the 
next power in the polar angle 9i. From simple scaling 
considerations it then follows that the corrections SH^ 
must be at least of order: 



5H2 
5H4 



0(a-i/2) 



(46) 



which all give marginal contributions to K2 for large as- 
pect ratio. 

The equilibrium value for a is obtained from the refer- 
ence free energy of the (undistorted) nematic. Combin- 
ing the ideal and excess parts of the nematic free energy, 
Eq. and Eq. ^ respectively, with the expression for 
Mo [Eq. (El])] gives: 

= (lnVp/(u) - 1) + ^ ((^;_,(ui,U2))) + ^. 

(47) 

On inserting the Gaussian averages for the ideal and ex- 
cluded volume contributions [Eq. (|^ ]. the asymptotic 
free energy becomes 



F 



GvdW 



ksTN 



\na + c^(7r/a)^/^. 



(48) 



where all terms independent of a have been omitted for 
compactness as they do not contribute to the degree of 
orientational order. Minimizing the free energy with re- 
spect to a yields the common quadratic form^^: 



(49) 



With this relation for a, an algebraic expression for the 
twist elastic constant for the SW spherocylinders can be 
formulated from Eq. ((36)) . After defining the spherocylin- 
der packing fraction as ~ {tt/4)xc (ignoring the end-cap 
corrections), we obtain 



K2D 
kBT' 



The twist elastic modulus primarily depends on the par- 
ticle packing fraction and the reduced temperature, with 
the aspect ratio merely playing the role of a linear scaling 
factor. 

Some remarks are now in order. The leading order 
contribution was found for infinitely thin roda^ and does 




0.5 



0.1 0.2 0.3 0.4 



FIG. 2; Twist elastic modulus K2 = K2Dx/kBT [Eq. (fSO]! ] 
for hard spherocylinders (l/T* = 0) as a function of packing 
fraction cj) = vop. The numerical ODF is given by the solution 
of Eq. p3|) . Inset: Temperature dependence of K2 for rods of 
aspect ratio L/D = 20 with attractive square- well (eooo > 0) 
and soft-repulsive square-shoulder (eooo < 0) interactions of 
range A = L -I- _D at a fixed packing fraction <j) = 0.2. 



not depend on temperature. The twist elastic modulus 
is independent of the range A of the SW potential. Cor- 
rection terms arising from the next leading order terms 
in Eq. ([^S]) are at least of order x ^ D/L and thus of 
marginal influence for sufficiently slender rods. Fig. [2] 
shows that the discrepancy with numerical results is very 
small even for relatively short rods with L/D = 10. The 
numerical data are based on a numerical evaluation of 
Eq. dH) using the exact ODF from Eq. It is im- 

portant to note that the cholesteric phase is only stable 
with respect to the isotropic state roughly when c ^ 1 or 
equivalently > x. As illustrated in Fig. [21 numerical 
results for the nematic solution of the ODF Eq. ([23)) are 
found only above a critical packing fraction. 

As expected, the twist elastic modulus is a monoton- 
ically increasing function of the packing fraction. The 
temperature dependence shown in Fig. [2] indicates an in- 
crease of the twist elastic resistance for attractive rods 
(^000 > 0), while the opposite trend is observed for a 
soft-repulsive square-shoulder potential (eooo < 0). The 
latter interaction may be particularly suitable for the de- 
scription of colloidal /d-virus rods as a crude model for 
the electric double layer or the polymer coat grafted onto 
the colloid surface^. 

The behavior of the twist elastic modulus with respect 
to the nematic order parameter is highlighted in Fig. [3] 
Since a - 0^ gg])] and a ~ 1/(1 - S") [Eq. ^] it 
is readily deduced that the asymptotic behavior of the 
twist elastic modulus for strong nematic order is given 
by the following scaling relation: 



KX 



1 



\-S 



1/2 



1 



3/2 



(51) 
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-0.45 -0.4 -0.35 -0.3 -0.25 -0.2 -0.15 -0.1 -0.05 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 

InS In1/(1-S) 

FIG. 3: Scaling of the twist elastic modulus K2 with respect to the nematic order parameter S. The results are obtained 
from the numerical ODF [Eq. (|23|l ] using L/D = 20. Continuous lines represent linear fits, (a) Quadratic behavior K2 ~ at 
moderate nematic order 0.6 < S < 0.9. (b) Asymptotic behavior K2 ~ 1/(1 — S)" with u > 0.5 in the regime of high nematic 
order (5" > 0.99). 



in terms of the constants a and b. This scaling relation 
is confirmed in Fig. [8)3 based on numerical results for 
the twist elastic modulus^. The crossover between the 
two scaling contributions in Eq. ([51]) is reflected by the 
fact that the exponent [1/(1 — S)'^] lies within the range 
0.5 < uj < 1.5 as 5 approaches unity (or equivalently 
1/(1 — S*)— >oo). A different scaling behavior is observed 
for moderate nematic order where K2 is found to be pro- 
portional to 5^ (Fig. [3^). This result is in agreement 
with previous theoretical predictions based on a Lcgcn- 
dre series expansion of the angular properties valid for 
weak nematic order— i^. 



rcos/3y + ^122 so that 

/ dri2 2;i2T2i2 = / dp (izi2Zl2 

Jx Jo J-x 

X / drrT2i2 
Jo 

= -A''(U1 • U2)(U1 X U2 • z). (53) 

The first spatial integral runs over the excluded volume 
of the spherocylinder for which we use the parametriza- 
tion advanced in Appendix A. This produces terms of 
0{LD^) and higher order in D which we will not show 
explicitly. We can thus write the first-moment integral 
as the following expression: 



B. Twist energy 

We proceed in a similar manner for the calculation of 
the twist energy. The first-moment spatial integral over 
the interaction kernel Eq. (PO]) is required first. None of 
the terms pertaining to the achiral parts of the interac- 
tion potential contribute to the twist energy so that 



Afi(ui,U2) = £212 / dri2Zi2T2i2 

-6212 j dYi2Zi2T2l2, (52) 



The second spatial integral is easily tackled by switching 
to cylindrical coordinates. Let us write ri2 = r sin /3x 4- 



A/i(ui,U2) = --e2i2A'^(ui • U2)(ui X U2 • i) + 0{LD^). 

(54) 

Inserting the asymptotic forms of for small polar an- 
gles into Eq. ([7]) and performing the azimuthal integra- 
tion over Lpi one can express the chiral torque as 




where the elementary Gaussian average {9"^ cos^ ip) ^ 1/a 
has been employed. This expression is similar to the one 
derived for infinitely thin rods^. The correction terms 
due to finite rod thickness are deemed to be very small, 
and can be ignored if the sphcrocylindcrs are not too 
short. 
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C. Cholesteric pitch 

A microscopic expression for the cholesteric pitch of 
SW sphcrocyhndcrs in the hmit of weak twist distor- 
tion is obtained by combining Eq. ([50]) and Eq. f55|) in 
Eq. (jl3p . After some rearranging one can write the pitch 
p = 27r/(7 in a convenient reduced form as 



where (+) refers to a square-well and (— ) to a square- 
shoulder potential. The factor e scaling the pitch com- 
bines the chirality parameter tc and the geometric pa- 
rameters of the range A and (inverse) aspect ratio x = 
D/L <^\. It can be interpreted as an integrated van der 
Waals energy which depends on the fourth power of the 
(SW) interaction range A: 

As with the twist clastic modulus, Eq. ([50]). the rescaled 
pitch depends only on the packing fraction and reduced 
temperature of the cholesteric phase. 

V. RESULTS AND DISCUSSION 

The density and temperature dependence of the 
cholesteric pitch obtained with the theory presented in 
the previous sections is depicted in Fig. [H The varia- 
tion of the pitch with packing fraction is non-monotonic, 
irrespective of the reduced temperature. The steep de- 
crease at low densities is a common result for slender 
rods. Since Kt oc (j)^ and K2 oc (j> [cf. Eq. ([55)1 and 
Eq. ([50)1 ]. the free-energy cost associated with the elas- 
tic deformation upon incrementing the density is more 
than offset by the simultaneous free energy gain due to 
an increase of the torque-field associated with the chi- 
ral interactions. The inverse proportionality p oa (j>~^ 
is in agreement with Odijk's result for rigid rods^ and 
consistent with experimental results of /d-rods at high 
ionic strengtbi^. At higher packing fractions, the orien- 
tational (nematic) order rises sharply and the influence of 
the finite rod thickness, embodied by the second term in 
Eq. (|56p . becomes apparent. The elastic resistance rises 
steeply and causes the cholesteric structure to unwind 
upon increasing density. 

For large temperatures (l/T* < 1), the behavior of 
the pitch is only weakly affected by the nature of the 
non-chiral interactions as can be inferred from the last 
term in Eq. ((56)) . From a qualitative point of view, the 
behavior is therefore the same for attractive square-well 
and soft repulsive square-shoulder potentials. Moreover, 
since p* cx T* , the temperature merely serves as a linear 



scaling factor and docs not influence the shape of the 
curve in Fig. 5^. 

It has been suggested that the increase of the pitch 
with density is due to pre-smectic fluctuations that coun- 
teract the twist deformationii^iSS. The present analysis 
would suggest that the observation can be accounted for 
within a simple mean-field theory for anisotropic yet ho- 
mogeneous systems in which one disregards density fluc- 
tuations; layered structures such as the smectic-A (SmA) 
phase possess inhomogencities in the average particle po- 
sition. That smectic fluctuations are not necessary to 
give rise to an increase in the pitch is also supported by 
experimental observations in colloidal /d-rods^^ where 
an unwinding of the cholesteric is observed at densities 
far below the cholesteric-smectic transition density. We 
stress that parallel configurations (induced by the finite 
rod thickness, i.e., spherocylindrical shape) that lead 
to an unwinding of the cholesteric state with increas- 
ing density also facilitate the formation of a smcctic-A 
phase. The two phenomena are therefore expected to 
be correlated. The possibility of a transition towards a 
smectic phase is not incorporated in our theory but it 
is anticipated that such an instability would lead to a 
much steeper increase and possibly a divergence of the 
pitch with density. Based on the phase diagram of pure 
hard sphcrocylinders^iiS^ the cholesteric-smectic transi- 
tion would occur at a packing fraction of (/) ~ 0.4, irre- 
spective of the aspect ratio (for sufficiently long rods). 
The unwinding of the cholesteric phase at high densi- 
ties is not unique to the fd system but has also ob- 
served in solutions of polysaccharide and polypeptide 
compound&i^ and DNAi^— . Theoretically, such a trend 
was first established in simulations based on the chi- 
ral hard spherocylindersS^, supplemented by a simple 
ground-state theory for a chiral hard Gaussian overlap 
model assuming perfect local nematic order—. 

The temperature dependence of the pitch is also highly 
non-monotonic. The near linear increase of the pitch with 
temperature for a system at constant density (or pres- 
sure) is a direct result of the achiral hard-core repulsion 
between the rods which dominates the (chiral) attrac- 
tive interactions at moderate to high temperatures. At 
T — > 00 the free energy of the system is governed en- 
tirely by the entropic contribution associated with hard 
spherocylinders, resulting in a nematic phase {p — > 00). 
The unwinding of the pitch with temperature has been 
found in solutions of polypeptidesi^ , and most notably 
for aqueous suspensions of /d- virus rodsi^ which was not 
accompanied by a thermal change in the intrinsic chiral- 
ity of the viral structure (e.g., due to a denaturation of 
the protein coating). Theoretically, this would translate 
to a chirality parameter ec which is constant, as we have 
assumed here. Our predictions of an increase in pitch 
with increasing temperature arc also consistent with the 
findings for the cholesterol ester CEEC— , a thermotropic 
mesogen which does not exhibit a smectic phase. 

The marked increase in pitch at low temperatures is 
in accordance with experimental observations in numer- 
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(b) 




T 



FIG. 4: (a) Scaled cholesteric pitch p* = {p/L)e [Eq. (|56p ] for attractive SW spherocylinders at reduced temperature T* = 
A;Br/|eooo| = 20 versus rod packing fraction (j> = vop- (b) Variation of the scaled pitch as a function of T* (on a logscale) 
for a constant packing fraction <j} — 0.35. Also shown is the behavior for a fixed rescaled pressure P* — Puo/|eooo| = 120 
corresponding to the isobar (thin continuous curve) in Fig. [T^. 




FIG. 5: (a) Phase diagram for attractive chiral SW spherocylinders {L/D = 10) of range X = L + D for various relative 
strengths of the chiral interaction = |£202/eooo|, where T* = kBT/\eooo\ is the reduced temperature and (j> — vop the packing 
fraction. CN denotes the cholesteric phase, (b) Temperature dependence of the reduced pitch p* = p/L of the cholesteric 
phase in coexistence with the isotropic. 



ous thermotropic systems the most ubiquitous being the 
derivatives of cholestero l^^'^° . A similar trend albeit less 
prominent is found for aqueous solutions of fd virus^^. 
The underlying scenario is analogous to the unwinding of 
the cholesteric with increasing density for the lyotropic 
case and originates from a stark increase in nematic or- 
der upon lowering T*. The twist elastic resistance as- 
sociated with the near-parallel configurations becomes 
anomalously large and gives rise to a strong unwinding of 
the cholesteric structure. As for the behavior with den- 
sity, an additional coupling and a divergent pitch is ex- 
pected at the nematic-smectic transition temperature as 
observed by Pindak et al^ and theoretically advanced 
by de Gennes^i. We do point out that in our particular 
model the decrease of the pitch with temperature (and 



constant density or pressure) observed at low tempera- 
tures is metastable with respect to an isotropic (gas)- 
cholcsteric phase separation (see the phase diagram in 
Fig.IIJi). 

Finally, from Eq. (j56p one can obtain a simple relation 
between the pitch and the particle aspect ratio. p/L oc x^, 
or equivalently p oc 1/L. The pitch of the cholesteric 
thus becomes tighter upon increasing the aspect ratio 
at a fixed mass density and temperature. This result is 
in line with observations in fd systems where the pitch 
is found to decrease for larger viral contour lengthsi^. 
The relation found in experiment, p (x l/L"'^^, reveals 
a much weaker dependency which could be attributed to 
the slight degree of flexibility, neglected in the present 
rigid-rod model. It is noteworthy that the Straley-Odijk 
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theor y'^^1 ° predicts the opposite trend p (X L. Their 
model is based on short-ranged (steric) chiral interactions 
induced by a thin hehcal thread of thickness D enveloping 
the rod. Here, the chiral interactions are long-ranged and 
scale with the rod length A oc i. A linear increase in the 
pitch with aspect ratio is predicted with a scaled Onsager 
theory for chiral hard Gaussian overlap rods with the ap- 
proximation of perfect local ncmatic order—. Although 
no microscopic justification for our scaling relation can 
be given, it does give the correct scaling of the pitch with 
contour length and is consistent with the long-ranged na- 
ture of the chiral forces between fd virus rodsi^. 

In Fig. [5] we shows that there is noticeable effect on 
the phase behavior of attractive SW spherocylinders for 
weakly chiral interactions. Upon increasing the strength 
of the chiral interaction relative to the dispersion forces 
(ec) the isotropic-cholesteric transition shifts to lower 
densities with an additional broadening of the two-phase 
coexistence region at moderate to high temperatures. 
Beyond a certain critical value, the isotropic gas-liquid 
envelope and the corresponding triple point become 
metastable relative to the direct coexistence between a 
low-density isotropic and a high density cholesteric phase. 
Since the chiral interactions are essentially attractive, we 
observe the same stabilization of the cholesteric state for 
repulsive square-shoulder potentials. The non-monotonic 
behavior of the pitch with temperature is reflected in 
Fig. [S] where the evolution of the pitch of the coexist- 
ing cholesteric phase is depicted. It is clear that pitches 
ranging from a few molecular lengths to a few hundreds 
of molecular lengths can be reproduced with our theory. 

The quantitative merits of the present theory can only 
be assessed by comparison with experiment for specific 
values of the interaction parameters of our model. The 
most important one is the chirality parameter Cc which 
captures the intrinsic strength of the chiral interactions. 
Apart from the shape of the mesogen, it is determined 
primarily by the nature of the long-range electrostatic in- 
teractions and the intricate surface structure of the par- 
ticles; in the case of viruses, for example, the chiral inter- 
action will depend on the helical configuration of surface 
charges. A change in temperature may induce confor- 
mational changes in the chiral structure which affect the 
magnitude or sign of the pitch. For example, a pitch in- 
version involving a sudden change of handedness of the 
cholesteric structure with temperature is known to occur 
in thermotropic systemsliS^. These issues are clearly be- 
yond the scope of our coarse-grained model and require 
a much more detailed representation of the molecular ar- 
chitecture of the mesoge n^^'^^ . 



VI. CONCLUSIONS 

An asymptotic analysis of the deformation free energy 
associated with a cholesteric phase of chiral spherocylin- 
ders is presented in this work. The rods consist of a hard 
spherocylindrical backbone with additional long-ranged 



achiral (attractive/soft repulsive) and chiral interactions, 
both represented by a simple square-well form with a 
range comparable to the length of the rod. Analytical 
expressions for the twist elastic constant and cholesteric 
pitch are deduced by invoking a Gaussian approxima- 
tion for the orientational distribution around the local 
nematic director. The approach is expected to provide 
an accurate representation of weakly twisted cholesteric 
states with a high degree of local nematic ordering which 
is essentially unaffected by the weak spatial variation of 
the director field. 

The results are relevant to both thermotropic meso- 
gens (e.g., derivatives of cholesterol) and more particu- 
larly lyotropic cholesteric systems such as fd virus rods 
where chiral interactions are mediated predominantly by 
long-range electrostatic forcesi^. The theory captures the 
behavior of the cholesteric pitch of fd rods, in particular 
its non-monotonic variation with temperature and den- 
sity, as well as the influence of the viral contour length 
(the effective aspect ratio). An extension of the th eory 
to simple hard-sphere chain models (e.g., see Ref. [971 ) 
would enable one to examine the effect of flexibility in 
more detail. 

For attractive SW rods, a steep increase of the pitch 
is found upon lowering the temperature, in line with ex- 
perimental observations in thermotropic systems. The 
steep unwinding of the pitch at low temperatures or high 
packing fractions is primarily due to a sharp increase of 
the local ncmatic order. The prevailing near-parallel rod 
configurations lead to an anomalous increase of the twist 
elastic resistance. This simple mean-field scenario con- 
trasts with the commonly expounded view in which the 
unwinding of the cholesteric is attributed to pre-smectic 
fluctuations which are geometrically incompatible with a 
helical structure. 

In the future, we plan to validate our theoretical flnd- 
ings with a simulation study of the current chiral sphe- 
rocylinder model along the lines of Refs. IosIIosIIoqI . This 
will allow us to test the accuracy of the Onsager-Parsons 
theory in predicting derivative properties such as the elas- 
tic constants of dense nematic systems. The simulations 
would also provide a better insight into the behavior of 
the cholesteric pitch close to a cholesteric-smectic transi- 
tion. Finally, it would be intriguing to study the impli- 
cations of the chiral interactions on the micro-structure 
of the smectic phase. 



Appendix A: Parametrization of the excluded 
volume of the spherocylinder 

The excluded volume manifold of two hard sphero- 
cylinders at flxed angle 7 is a spheroparallelepiped (see 
Fig. [5]) which is most conveniently parametrized by 
switching from the laboratory frame to a particle frame 
based on the orientational unit vectors Ui . Let us further 
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define the unit vectors 



Ui X U2 

I sin7| 
Uj X V, (i=:l,2) 



(58) 



so tliat {ui, V, Wi} are two ortlionormal basis sets in 3D. 
The centre-of-mass distance vector can be uniquely de- 
composed in terms of these basis vectors: 

ri2 = (ri2 ■ Uj)u,, + (ri2 • v)v + (ri2 • Wj)wi. {i = 1, 2) 

(59) 

The leading order contribution to the excluded-volume 
body is of 0{L^D) and stems from the overlap of the 
cylindrical parts of the spherocylinders. This gives rise 
to a 3D parallelepiped (central section in Fig. [6]) which 
can be parametrized as 



L 



-tiUi 



-t2U2 + -DtsV, 



(60) 



with —1 < ti < 1 for I = 1,2,3. The Jacobian as- 
sociated with the coordinate transformation is Jcc = 
■jL^Dj sin7|. Higher order contributions account for end- 
cap effects due to the finite thickness of the rods^^S. The 
first correction term is of 0{LD^) and originates from an 
overlap of the hemispherical end-cap of one rod with the 
cylindrical part of the other. The resulting half cylinders 
on the four edges of the parallelepiped, at the boundaries 
of the excluded volume in Fig. [SJ can be parametrized as 
follows: 



"^12 

^HC± 
'12 



=F^U2 + "^^i^i + i2-D(±wi cos/3 + vsin/3) 

±-^Ui + ■^iiU2 + t2D{±W2 COS/3 + vsin/3), 

(61) 



with boundaries — l<ti<l,0<t2<l and tt/2 < 
(i < 7r/2 and Jacobian Jch(hc)± = \LD'^t2- The final 
contribution stems from the spherical segments located 
at the four corners of the parallelepiped (yellow sections 
in Fig. [5]). These arise from the overlap of two hemi- 
spherical end-caps and are of 0{D^). The segments can 
be parametrized by invoking a different orthonormal set 
{u+, u_, v} where 



U± = Ui ± U2, U± = 



|u±| 



(62) 



which leads to the following expression for the centre-of- 
mass distance corresponding to the four segments: 



'l2 



'12 



± — u+ + rD (cos + sin 9 sin </'+ fl- 
it sin 6' cos </3_u_(_) 

±— u_ + rD{cos Ov + sin 9 sin ip-U+ 



± sin 6* cos (y9_|-U_), 



(63) 



with parameter intervals 0<r< 1, < 9 < tt, (7 — 
7r)/2 < ip+ < {n - 7)/2 and -7/2 < (p- < 7/2. The 




FIG. 6: (Color online). Illustration of the excluded vol- 
ume manifold of two spherocylinders for a relative orien- 
tation angle 7 and the unit vectors defined in Eq. (|58|l . 
The cylinder-cylinder CC, cylinder-hemisphere CH, and 
hemisphere-hemisphere HH volumetric sections are indicated 
in green, red and yellow, respectively. 



Jacobian associated with the coordinate transformation 



is J, 



{HH+±){HH-±) 



D^r"^ sin( 



For the twist clastic modulus we require the second- 
moment excluded volume [cj. Eq. (j29p ]: 



(64) 



?;m2(ui,U2) / (iri2(z • ri2)^ 



which can be evaluated separately for each of the ex- 
cluded volume sections using the parametrization above. 
The integrations over the parametrization variables can 
be worked out without difficulty. Let us define the fol- 
lowing dot products 





= {L/2)xi, 


B 


= Dwi 




= -Dwi • z 


E± 


= Dvi± ■ z 



(65) 



The cylinder-cylinder CC contribution then reads 



ii^£)|sin7| J dti J dt2 J dt3{z ■ rccY 



'-L^D\su-i-f\{Al + A^ + B^). 



(66) 
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For the cylinder-hemisphere CH contributions we obtain: 



7r/2 



d(3 



1 J<3 J-it/2 

{{i ■ + ■ vchS + ■ ■ ■] 

LDH^-^{A,C,^A,C,) + ^{Al+Ai) 



(67) 



Finally, for the hemisphere-hemisphere HH contribu- 
tions we have: 



the subsequent angular averaging. It is enlightening to 
expand the argument as a Taylor series: 



2tl 



■a; + 0(x3), (72) 



4x2 + 4^2 +13 + 

from which one finds that the leading order CC contri- 
bution is of 0{LD^) and thus of marginal importance 
for the relevant range of aspect ratios L/D > 10. The 
integration over the CH parts requires more effort, but 
the result can be cast in a similar form: 



= iiD25|sin7|cos7J-c//(a:,7), (73) 



/ drr^ / dOsine 
'0 Jo 

n(ir-7)/2 

dip+{{z-rHH+^)^ + ■■■} 

/(7-7r)/2 



/ dip^{{z-rHH^J' + ■■■} 

J-j/2 



= DUll{E.\n.\r + ^{7r~j){E+\u+\r 



TT I |u+|i?:j: cos 



7 



|u_|£:^sin? 



}. (68) 



with 

, ,^ sin^ /3 

^CH[X,J) = I dti I dt2t2 / dfj — CH^- 



tt/2 
tt/2 
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Here, r'^^^ represents the vector norms: 



(74) 



^f2"" = V4^ + ^+*i±|2+|cos/3|sin7|. 

(75) 

Without further analyzing Tch it is evident that the CH 
contributions are small and at least of 0{LD^). The HH 
integrations will produce terms of even higher order in 
D and therefore the analysis need not be pursued any 
further. 



It is easily verified that all contributions are symmetric 
under inversion Ui — Ui and interchanging Ui ^ U2, 
as required. The total expression Eq. ([64|) is obtained 
by adding the contributions from the different sections, 
-I- , which provides an analytic 



VAh 



CC 

-'a/2 



result for the moment excluded volume of two sphero- 
cylinders of arbitrary aspect ratio. 

For the twist energy one requires the first-moment ex- 
cluded volume over the pseudo-scalar T212. Recalling 
Eq. (dnD we get 

?;mi(ui,U2) = C0S7 / dri2(z-ri2)(ui xu2-ri2), (69) 

with cos 7 = Ui • U2 . Exploiting the orthogonality of the 
unit vectors and collecting terms gives for the leading 
order contribution: 



1 



"Ml -^L-DBsin cos Tec (xn)- (70) 



Here, J-cc represents a triple integral: 



3 »i 

-^cc(a;,7) = n / 
i=i -^"1 



dt, 



^3 



+ ^+t% + ^f^cos7 



(71) 

which cannot be solved in closed form. In the asymptotic 
limit (cos 7 '--^ 1) Tec becomes a function of the inverse 
aspect ratio x = D/ L only, and no longer plays a role in 



Appendix B: Gaussian averages 

The Gaussian averages required for the evaluation of 
Hq and H2 [Eq. and Eq. (gl])] have been deduced by 
Odijk^^. We quote them here: 



\i\om+oi))) 



^^^1/2^ -5/2 

2 

21^1/2^-5/2 

2 



(76) 



The last term in Eq. (j44|) can be evaluated in two steps. 
Since 6*^(^1 — 6'^)/|7| ~ 0{9'^) a simple Gaussian integral 
suffices to establish the following scaling result: 



I7l 



(x a I dOO exp 





pa 



-3/2 
. 

The pre-factor p is obtained by numerical evaluation of 
the following triple integral: 

(.ir/2 i'-k/2 p2iT 

p = / ^6*1611 / d9202 



dA(p 
'2^ 



a 



I7l 



X exp 

which gives p ~ —1.772 for any given a ^ 1. 



(78) 
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